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SUMMARY 


The characteristics of the position ’and form of 
the transition surface through the critical velocity 
are computed for flow through flat and round nozzles 
from subsonic to supersonic velocity* Corresponding con- 
siderations were carried out for the flow about profiles 
in the vicinity of sonic velocity. 


I. INTRODUCTION 


While useful methods of calculation exist for pure 
subsonic and for pure supersonic flows of compressible 
gases difficulties 'arise In the mathematical treatment 
of mixed flows with subsonic and supersonic ranges. Such 
mixed flows exist; (l) in the transition from subsonic to 
supersonic speeds in Laval nozzles, (2) about profiles 
in a flow at high subsonic speed on the appearance of 
local supersonic ranges, (3) in front of blunt bodies 
in a flow at supersonic velocity with a local subsonic 
region behind the compression shock in the vicinity of* 
the stagnation point. 

The present report takes up the first as the 
simplest of the three problems named. ' It is a fact that 
for this case a rough view of flow conditions suffi- 
cient for many problems can be obtained by the methods 
of hydraulics, that is, the nozzle is considered a 
flow tube and the magnitude of the velocity Is considered 
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constant in every cross-sectional plane. The critical 
velocity (flow velocity = sonic velocity) is then reached 
exactly in the minimum cross section of the Laval 
nozzle. In the strict two-dimensional treatment, 
however, a curved line (critical curve) is obtained for 
the passage of the velocity through the critical value 
in the plane of the longitudinal section of the nozzle; 
it begins at the nozzle wall in front of the minimum 
cross section and cuts the nozzle axis downstream of 
the minimum cross section. (See fig. 1.) 

In various reports, especially those of Th. Meyer (1) 
G. J. Taylor (2), H. Gortler (J), and Kl. Oswatitsch and 
K. Hothstein (II), expansions in power series for the 
determination of the critical curves are given. In 
what follows, general statements on the position and 
curvature of the critical curve for « given nozzle are 
obtained from such expansions m power series by 
breaking off the series after the first two terms. The 
investigations are concerned with flat and round nozzles 
and provide sufficiently accurate results for practical* 
pvirposes for nozzle curvatures that are not too sharp. 

Corresponding expansions in power series are 
applied, in conclusion, to the flows through flat 
nozzles with curved axes. In this way, information is 
obtained on the variation of the critical curve for 
profile flows with local supersonic regions. 


II. POTENTIAL EQUATION 


The nozzle axis is selected as 
origin of the coordinate system Is 
on the axis at which the critical 
reached, (See fig. 1.) Tith the 
vortical flow and perfect flow of a 
constant specific heats c , c v (k 


the x-axis and the 
placed at the point 
velocity Is first 
hypothesis of non- 
perfect gas with 
“ c /c ) the potential 

ir v 


equation 


du 

Ox 



- ) + 




2uv 




0 


( 1 )' 
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holds. In this u, v are the x- and y-components of the 
flow velocity and a, the local sonic velocity, which 
is related to u, v and the critical velocity a ,* through 


p k + 1 iO k - 1 

j — a** - — 2 — V u + v y 


o=0 characterises the flat or two-dimensional flow, 
0=1 the three-dimensional flow through Btresm tube3 -of 
of circular cross section; x, y are cartesian for O = -0, 
cylindrical coordinates for a = 1. 

Substituting (2) in equation (1) and introducing 
the dimensionless velocity components 


U = u/a*. 


V = v/a* 


the potential eq\.xation (1) becomes; 


(, .,2 k - X „-}\ .6 V /, k - 1 „? „?\ 

S ^ J ' FTT ‘ J + 3? 1^ " kTT u ' v J 


— k UV + a l - tt-— 

k + 1 dy uv 1 lc + 1 


(b 2 + V 2 ) = 0 (1+) 


Since the present investigation was limited to the 

vicinity of the critical curve, 

» 

U = 1 + u, V = v (5.) 

in which u and v are small quantities, and equation (iq) 
becomes 


du . ? , k - 

“ \2u + u^ + iT -i 


J] _ 

+ 1 v J c5y k + 1 


(k - 1) 


?• u 2 - v 2 + t — Lp (1 + u) v - a r — —■ 

lc + 1 k + ,1 Oy k + 


* - iH4 r + A f=°- 


k + 1 
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As x—^O, y — >0, u, and v approach zero and also 
on the basis of (6) 6v/dy and, consequently, the 
quotient v/y approach zero, while the velocity 
rise 6u/dx along the axis ought not vanish at the 
origin. Then if the small quantities u, v, 6v/6y, 
and v/y are considered linear only, the following 
approximate relation is obtained from (6) 


! (k + 1) u 




+ 2v 


ou 

r -'y 



( 7 ) 


Since the nozzle flow is symmetrical with respect 
to the x-axis, ou/8y also approaches zero for x~>0, 
y — . Consequently .the term 2v 6u/dy may be ignored, 
by means of which ( 7 ) becomes further simplified to 


i /, , , s ou ov _ v 

{ k + 11 u 5J " ’ 0 y 

f 



( 8 ) 


*'/ith o=0, (8) was already applied in another con- 

nection by Kl. Oswatitsch and Ii. Oieghardt ( 5 )* 


III. FLAT NOZZLES 


Equation (8) furnishes general relations, for a = 0, 
for the flow through flat nozzles in the vicinity of the 
critical curve. Like Kl. Oswatitsch and W. Rothstein ([(.), 
taking into account the symmetry around the x~axis, 
setting 


= f Q (x) - 1 - y 2 f 2 (x) + y^+f^(x) +... 


u 


= = f’ 0 (x) + y 2 f',(x) + y^-f ’ |, (x) ... 

^ ■ X. 




( 9 ) 


V = = 2yf (x) + iiy3f. (x) + 


8y 


4 


* 1 1 


J 



NACA TM No. 11^7 


5 


in ivhich the primes signify derivatives with respect 
to x, and, obtain from (8). 


(k + 2) (f + y 2 f' + ...) (f» + y 2 f" +... ) 

0 2.0 2 

= 2f 2 + 12y 2 f^ +. . . + 2 a (f 2 + 2y 2 f^ + . . . ) 

By arranging in order of powers of y and equating 
coefficients of the individual powers to zero this 
system of equations is obtained 


2(1 + a) f = (k + 1) f» f" 

2 0 0 

2(6 + 2a*) ft = (k + 1) (f» f« + f" f» 

4 v 0 2 0 2/ 


( 10 ) 


through which the functions f (x), fj (x)... in order 

cl * 4 - 

can be expressed by the function f' (x) and its 
derivatives. The function 

f 'o = u o (x) 

which characterizes the velocity distribution along the 

nozzle axis, remains undetermined. Given u (?) the 

’ o 

flow is completely established through (10) and can be 
understood to be nozzle flow, if a pair of streamlines 
that are mirror images of one another with respect to 
the x-axis are assumed as nozzle walls • For further 
discussion, setting 


f 'o = u o (x) 


a x 


( 11 ) 


from which 


f = 
2 


k + 1 

2 


a 2 x. 



(k + l) 2 
2 % 



(12(a)) 


"Translator's note: The number 2 was omitted from the 

German report. The correction here changes several 
subsequent equations and several values in the table. 
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follows Immediately from (10 ) with O = 0 . According 
to (9) 


u = ax + a 2 y 2 + 


v = (k + 1) a 2 xy + — a 5 y 3 + , , . 


is obtained for the components of the flow velocity. 

'xhe dots represent terms of higher order of y. If ( 11 ) 
is taken as the first term of a power series expansion 
for u Q (x), the expressions for f 2 and f^ in (12(a)) 

are also the first terms of power series expansions. 

From ( 13 (a)) the critical curve Is obtained from the 
requirement 

• ' • 1 ‘ ‘ ' ’ ? 2 

(1 + u) + v = 1 

from which in (13(a)) the parabola 

x K = - “y K 2 (lMa)> 

is obtained, therefore, with u =0 as a first 
approximation. It. cuts the nozzle axis, normally, at 
the point x = y = 0 and its curvature there is 

l/p K = (k + 1) a (15(a)) 

The vertices, therefore the points, of the streamlines, 
adjacent to the nozzle axis, are given by v = 0 and, 
according to (13(a)), lie on the curve 

Xo = - k t' -- ay 2 (16(a)) 

•-> O O 

Considering, now, the streamlines at a given vertex 
distance y^ (fig. 1) from the nozzle wall, then 
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is the distance from the center .point of the narrowest 
cross section downstream to the intersection point of 
the nozzle axis and the critical curve. 

The junction point of the critical curve with the 
nozzle wall is upstream and is obtained, as a first 
approximation from (li^a)) with y^ = yg, which yields 


11 = -(*K - *s) = ^ ay s 2 = ££ 

for the distance rj from the junction point to the 
narrowest cross section. 


Finally the calculation of the curvature at the vertex 
of the nozzle wall is to be made. In moving outward 
from the vertex to the wall by an element of 
curve ds = dx, the tangent turns through the angle 


dm 


d v 

1 + u 


1 dv 
1 + u 3 x 


dx 


1 ov 

r+u sx 


ds 


Therefore the curvature • at the vertex is 

1 dm I 6v 

pg dS “ 1 + U ?x 

and taking into account (13(a)) and (15(a)) 

l/pg = (k + 1) a 2 y s = l/p K ay s (19(a)) 

Up to now the velocity distribution u D (x) had been 
considered as given by equation (11) and the nozzle 
walls computed for it. Practically,' the reverse is 
the problem, namely to ascertain for a particular 
nozzle, that is for assigned values of y s , pg, the 
flow in the narrowest cross section, therefore the 
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values a, rj, and p, * On the basis of the equations 
that have preceded it xs possible to write down the 
solution to the problem, immediately, namely 



(20(a)) 


If the flow Is in the positive direction of the x-axis, 
as In figure 1, then the subsonic region is to the left 
and the supersonic region to the right and, accord- 
ingly, a> 0. 

In practice it is recommended that all distances be 
expressed, dimensionless, in terms of y^, and 

accordingly set y_ = 1 throughout (20(a)). 

O 

IV. ROUND NOZZLES 


Corresponding relations for round nozzles come 
from ($) and (10) with 0=1, Retaining (11) for 
the velocity distribution u q (x) along the nozzle 

axis, the following relations are obtained in place 
of (12(a)) through (20(a)): 



k + 1 p 



(k + 1) 

£5 



(12(b)) 


NACA TM No. 1147 


9 


k + 1 p p 

u = ax + — j- — a^y- + . . 


k + 1 ? 

v = — 2 — orxy 


(k + 1 'T?,*. 


TF~ ^ + -** i 


(13 Co)) 


k + 1 p 

X K = ~ ~ a ^ r K 


(li).(b) ) 


1 k + 1 


P K 


a 


(15(b) ) 


k + 1 ? 

x s = 8~ aj s 


( 16 (b)) 


k + 1 p 
e = -x 0 = — - — ayo^ 

o 8 ° 


(17(b) ) 


il = -(x K - x s ) = =■ (k + 1) a y s 2 = e (18(b)) 


1 _ k + 1 _2 


a* -y s 




(19) 
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a = 


*K = 


i 2 

j (It + 1 ) p g y s 


J 2p q y^ 
i ic + r 


■i] = e = \j 2 (k + 1) 


8 


(2C(b) ) 


To compare a flat nozzle with a round, one of the same 
longitudinal section, formulas (20(a)) and (20(b)) 
are compared. (See fig. 2.) 


For the flat 

nozzle' rj = 

2€, f 

Moreover 

6 round 
e flat 

* i 

r* V 2 

4 ■ 


T) round 

5 ,r; 


"0 flat 

o \ - 


a round 
a flat 

y 2 = 

The velocity 

increase is. 

there: 


nd 


~ A r’ ~Z 

0 . op 


i±l 


larger for the round nozzle than for tuc flat nozzle. 
All of these results are fully confirms a by the nozzle 
flows calculated by Rl. Oswatitsch and Ao that sin, as 
the enclosed table on page i3 shows. 
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V. COMPARISON / m TIJT PLOW- tub:: .iPPhOXIJTiflON 


The author's results are r. ow to be comer-. re l with 
the flow- tube approximation by CompuLink the increase 
in velocity a by formulas (20(a)) and (20(b)), and 
by the flow- tube cheery. For this purpose the flow 
density along the nos sic axis in the vicinity of the 
narrowest cross section is expanded, as follows 


8 = XT 
P* 

N _L._ 

/ lr + 1 

= const. U h . . 

\1C - J. j 


■ 1 __ 
«2\lc-1 


~ consc 

~ const. (1 - u”) 

Taking the velocity dis tribu cion (13) as a bo sis, for 
the nozzle walls in the .neigh' -orb . j .d of the narrowest 
cross section the following is obtain-:. 3. in the case of 
a flat nozzle 

-1 

/ _ s' \ f - '2 1 

77 


. / V 

+ 1 

- 1 - 


u) \r 

T 

2u 




3 

u) !r 

2 

- 1 

- 2u) ] 
/ 

:-l 

u 2 ) 




(l - u 2 ) = const. = y^, it ~ j (l - u 2 ) = (l + a 2 x 2 ) 

and in the case oF a round noszle 


Fh - 


b 


u l ~ consc. 

j 


2 


S 


> j 


1 

( -A"? 

y (1 - U g = 


A 1 + 


o o\ 
a 2 '- 1 
~ j - < 


From this , the curvature at the vortex is 

p 


and 


i/p q = 


i/p„ = a2 s, 


respectively , 
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and in place of the first equations of (20(a)) 
and (20(b)) 


a 


V 2 p s y s 


(2.1(a) ) 


a 



(2.1(b)) 


The comparison of (20(a) ,20 (b) ) and (21 (a ) ,21 (b) ) shows 
that the velocity increase a in the flow-tube approxi- 
mation fo r both flat and round nozzles is too large by 

a factor ^ — ~ , that is for k = 1 • i , about 10 percent. 

The values for the flow-tube approximation are added 
to the table, page l8, in square brackets. 


VI. FLO ,7 ABOUT PROFILES 


Through similar power-series expansions the flow 
past a profile in the vicinity of transit through the 
critical velocity can be discussed. If a streamline 
adjacent to the profile is thought of a rigid wall, 
the flow may also be explained as flow through a 
nozzle with a curved axis. The origin of the coordinate 
system is located at that point of the profile at 
which the critical velocity is just reached, the 
x-axis downstream on the profile tangent and the y-axis 
normal to the profile, outward, (See fig. 3 ») 


Since the flow is no longer symmetrical to the x-axis, 
instead of (8) the rather complicated equation (7) 
which must be specified with 0=0, by all means, is 
taken as a basis and in place of (9) 


cp = f Q (x) + yf^x) + y^f .p (x) + y?f (x) +... 
u = || = f» 0 (x) + yf 1 x (x ) + y 2 f ' (x) + y 5 f' 5 (x) 

V = = f (x) + 2yf (x) + 3y 2 f (?) +.,. 

o y 1 2 . 3 


+ 


• • • 


V 


J 


( 22 ) 
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By- putting (22) in ( 7 ) 

(k + l)(f» Q + yf (f” 0 + yf"^..*) = (2f 2 + Syf +...) 


and by comparing coefficients 


f 2 = f 0 f" 0 + f x 


of-. 

2 




= (lr + 1 ) (f f" + f 11 f 1 ^ + k/f f 1 + ft r 

■ c L) o 1 1 x cr 1 ) 2 i z) 

J 


( 23 ) 


while only the function f’gCx) had remained arbitrary 
in the nozzle flow, here there are two arbitrary 
functions f’ 0 (x), f-^x). dith f' 0 (>') and f 1 (x) 

the remaining coef f ic iont func fcior.s f £ (x ) , f >;(>:),.. . 
in order may be computed from the system of eyuati ons ( 23 ), 
by which the flow is completely established. 


Similar to (11) set 


0 

and obtain from (25) 


f» = ax, fp = -px 


(2 k) 


^Z~ a ^’ + P 2 ) x 


r 3 = - S_+_i a|3 - a 2 + P 2 )x 


and from (22) 
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u = ax - (3y + 

1"— a 2 + p 2 ) y 2 +... 

v = - (3x + 

(k + l)a 2 + 2(3 2 ; 

1 

— 4 

—1 

xy 

j 

- < iL -+-la(3 + 14.(3 

(k + 1) „Z A rj 2 

— a + p 

1 2 

x x* y + • • « 
} 
i 

1 


Here a signifies the velocity increase along the profile 

on transit through the critical velocity; it is therefore 
positive at the starting point A and negative at the 
terminal point 3 (fig. 4) of a local supersonic region. 

The profile curvature is given with (3 > 0 at the 
point A or B as 





(26) 


For the slope of the critical curve 'with respect to 
the profile tangent with 

(1 + u) 2 + v 2 ~ 2u “ 0 

from ( 25 ) the relation 


tan & 




(27) 


is obtained. 

'"'Translator’s Note: This 2 was omitted through 

error in the original German report. The correction 
applied here changes values of the constant in the 
succeeding equations and some values in the table. 
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For tha slope of the streamline at a distance y., 
from the profile 0 



s P - j^(k + 1) a 2 + j3 2 J y g 
is obtained -and taking ( 26 ) into consideration 



(23) 


according to (27) the critical curve at A (a > 0) 
and B (a < 0) is steeper the larger the velocity 
increase or decrease and the less the profile is 
curved. By (2o) in the limiting case of a = 0 the 
circle of curvature st S is concentric with the 
circle of curvature of the profile and, as is already 
known from- the- flow-tube approximation, becomes still 
flatter with increasing a. 


VII . SUMOhY 


Approximation formulas were developed for the 
position and curvature of the critical curve for 
trans ition through the critical velocity in the 
neighborhood of the narrowest cross section of flat 
and round Laval nozzles. In comparison with 
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nozzle flows calculated by Kl, Oswatitsch and 
W* Roths te in (J+) they showed satisfactory agreement* 

In addition, corresponding approximation formulas were 
deduced for the flow at profiles with local super- 
sonic regions. 


Translated by Dave Peingold 
National Advisory Committee 
for Aeronautics 
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COMPARISON OP APPROXIMATION FORMULAS (20(a)) AND (20(10) 
WITH THE NOZZLE PLOWS COMPUTED BY KL. OSWATITSCH 

AND W. ROTHS TEX' 'I ' 


10 


p s /y s = 5 


10/3 


a flat ! 

i 

[a. f lo.tTj I 
! r Tube J j • 
a "round, j 

a rouncQ i 
_ Tubs j i 


e flat 
e round 


T) flat 
T) round 


0 .20 ( 0 . 20 ) 
[ 0 . 22 ] 

0.29 (0.28) 
[°- 52 i 

0.2J (0.27) 

[<>• >3 

. O.It-l (0.37) 

' r 1 
j 0 > 4 pj 

• 0 . 35 . ( 0 . 32 ) 

ip. 59 ] 

0.30 (0.I4.2) 

- [0.55] 

0.08 (0.08) 

0.12 (0.12) 

o.ii. (0.16) 

0.09 (0.10) 

0.12 (0.li+) 

O.ll; ( 0 . 13 ) 

0.16 (0.16) 

0.25 ( 0 . 25 ) 

0.20 (0.25) 

0.09 (0.08) 

0.12 (0.12) 

0.ll|. (0.14) 


The unbracketed numbers 'have been computed y/ith approxi- 
mation formulas (20(a)) and (20(b)), the numbers in 
curved brackets taken from figures 7 and 10 of the. 
paper by Kl. Oswatitsch and //, Rothsteln (4). The 
square bracketed numbers are from formulas (21(a)) 
and (21(b)) for the flow- tube approximation. 
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Figs. 1,2 


P 



Figure 1. Illustration of terms 
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Figs. 3,4 
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Figure 4. Profile with supersonic region. 


